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Abstract. The Hamiltonian H specifies the energy levels and the time evolution of 
a quantum theory. It is an axiom of quantum mechanics that H be Hermitian because 
Hermiticity guarantees that the energy spectrum is real and that the time evolution 
is unitary (probability preserving). This paper investigates an alternative way to 
construct quantum theories in which the conventional requirement of Hermiticity 
(combined transpose and complex conjugate) is replaced by the more physically 
transparent condition of space-time reflection {VT) symmetry. It is shown that if 
the VT symmetry of a Hamiltonian H is not broken, then the spectrum of H is 
real. Examples of T'T-symmetric non-Hermitian quantum-mechanical Hamiltonians 
are H = + ix^ and H = — x^. The crucial question is whether T'T-symmetric 
Hamiltonians specify physically acceptable quantum theories in which the norms of 
states are positive and the time evolution is unitary. The answer is that a Hamiltonian 
that has an unbroken VT symmetry also possesses a physical symmetry represented by 
a linear operator called C. Using C it is shown how to construct an inner product whose 
associated norm is positive definite. The result is a new class of fully consistent complex 
quantum theories. Observables are defined, probabilities are positive, and the dynamics 
is governed by unitary time evolution. After a review of T'T-symmetric quantum 
mechanics, new results are presented here in which the C operator is calculated 
perturbatively in quantum mechanical theories having several degrees of freedom. 
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1. Introduction 

In this paper we present a brief review of some recent work on an alternative way 
to formulate of quantum mechanical models and present new results concerning the 
perturbative calculation of what has become known in this theory as the C operator. 

In the conventional formulation of quantum mechanics the Hamiltonian H, which 
incorporates the symmetries and specifies the dynamics of a quantum theory, must 
be Hermitian: H = W . The usual meaning of the symbol f, which indicates Dirac 
Hermitian conjugation, is combined transpose and complex conjugation. It is commonly 
thought that a Hamiltonian must be Hermitian in order to ensure that the energy 
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spectrum (the eigenvalues of H) is real and that the time evolution of the theory is 
unitary (probability is conserved in time). Although H = W is sufficient to guarantee 
these properties, it is not necessary. Indeed, we believe that this condition of Hermiticity 
is a mathematical requirement whose physical basis is somewhat obscure. Recently, a 
more physical alternative axiom called space-time reflection symmetry {VT symmetry), 
H = H^^ , has been investigated. This symmetry allows for the possibility of complex 
non-Hermitian Hamiltonians but still leads to a consistent theory of quantum mechanics. 

Because VT symmetry is an alternative condition to conventional Hermiticity, it 
is now possible to construct infinitely many new Hamiltonians that would have been 
rejected in the past because they are not Hermitian in the usual sense. One example of 
such a Hamiltonian is H = ^p^ + ^fi^x^ + iex^, which is the quantum mechanical analog 
of the quantum field theoretic Hamiltonian H = J die [^tt^ + |(Vv2)^ + |/^^V^^ + ie^^] ■ 
Another example of a PT-symmetric non-Hermitian Hamiltonian is H = -|- i/i^x^ — 
ex^, which is the PT-symmetric analog of the quantum field theoretic Hamiltonian H = 
J (ix [ivr^ + i(V(y5)^ + l/i^v'^ — This latter Hamiltonian could be an interesting 

candidate for describing the Higgs sector of the standard model. It should be emphasized 
that we do not regard the condition of conventional Hermiticity as wrong. Rather, we 
view the condition of VT symmetry as offering the possibility of studying new kinds 
of quantum theories that have heretofore never been studied because they have been 
thought to be physically unacceptable. 

Let us review the properties of the space reflection (parity) operator V and the 
time-reflection operator T: V is a linear operator with the property that V^ = 1 and 
has the effect p —p and x —x; T is an antilinear operator with the property 
that = 1 and has the effect p —>■ —p, x ^ x, and i —i. The operator T is 
called antilinear because it changes the sign of i. We know that it reverses the sign 
of i because, like V, this operator preserves the fundamental commutation relation of 
quantum mechanics, [x,p] = i, known as the Heisenberg algebra. 

It is easy to construct Hamiltonians of the form H = p'^ + V{x) that are not 
Hermitian but do possess VT symmetry. The trick is to take the potential to be a 
function of ix: V = V{ix). We also impose a general condition that has not been widely 
emphasized in the literature; namely, we require that H be symmetric: H = H^, where 
T represents the transpose. The reason for this symmetry condition will become clear 
later on. For example, consider the one-parameter family of symmetric Hamiltonians 

H = p^ + x^{ix)\ (1) 

where e is real. While if in (0) is not symmetric under V oi T separately, it is invariant 
under their combined operation. Such Hamiltonians are said to possess space-time 
reflection symmetry. Other examples of complex Hamiltonians having VT symmetry 
are H = p^ + x^iixY, H = p^ + x^iixY, and so on pp. Note that these classes of 
Hamiltonians are all different. For example, the Hamiltonian obtained by continuing H 
in along the path e : 8 has a different spectrum from the Hamiltonian that is 
obtained by continuing H = p^ + x^{ixY along the path e : — > 4. This is because the 
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Figure 1. Energy levels of the Hamiltonian H ~ + x'^(ixy as a function of the 
parameter e. There are three regions: When e > 0, the spectrum is real and positive 
and the energy levels rise with increasing e. The lower bound of this region, e = 0, 
corresponds to the harmonic oscillator, whose energy levels are En ~ 2n + 1. When 
— 1 < e < 0, there are a finite number of real positive eigenvalues and an infinite 
number of complex conjugate pairs of eigenvalues. As e decreases from to —1, the 
number of real eigenvalues decreases; when e < —0.57793, the only real eigenvalue 
is the ground-state energy. As e approaches — 1+, the ground-state energy becomes 
infinite. When e < — 1 there are no real eigenvalues. 

boundary conditions on the eigenfunctions are different. 

Tlie class of PT-symmetric Hamiltonians is larger than and includes real symmetric 
Hermitians because any real symmetric Hamiltonian is "PT-symmetric. For example, 
consider the real symmetric Hamiltonian H = p'^ + + 2x. This Hamiltonian is time- 
reversal symmetric, but according to the usual definition of space reflection for which 
X — > —X, this Hamiltonian does not appear to have VT symmetry. However, the parity 
operator is defined only up to unitary equivalence, and if we express the Hamiltonian 
in the form H = p"^ + {x + 1)"^ — 1, then it is evident that H is VT symmetric provided 
that the parity operator performs a space reflection about the point x = —1 rather than 
X = 0. See Ref. |2] for the general construction of the relevant parity operator. 

With properly defined boundary conditions the spectrum of the Hamiltonian H in 
is real and positive when e > |3] and the spectrum is partly real and partly complex 
when e < 0. The eigenvalues have been computed numerically to very high precision, 
and the real eigenvalues are plotted as functions of e in Fig. [TJ 

The VT symmetry of a Hamiltonian H is said to be unbroken if all of the 
eigenfunctions of H are simultaneously eigenfunctions of VT. Note that even if a 
system is defined by an equation that possesses a discrete symmetry, the solution to 
this equation need not exhibit that symmetry. For example, although the differential 
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equation y(t) = y(t) is symmetric under time reversal t —t, the solutions y(t) = e* and 
y{t) = e~* do not exhibit time-reversal symmetry; other solutions, such as y{t) = cosh(t), 
are time- reversal symmetric. The same is true of a system whose Hamiltonian is VT 
symmetric. Even if the Schrodinger equation and corresponding boundary conditions 
are VT symmetric, the wave function that solves the Schrodinger equation boundary 
value problem need not be symmetric under space-time reflection. When the solution 
exhibits VT symmetry, we say that the VT symmetry is unbroken. Conversely, if the 
solution does not possess VT symmetry, we say that the VT symmetry is broken. 

It is extremely easy to prove that if the VT symmetry of a Hamiltonian H is 
unbroken, then the spectrum of H is real: Assume that (i) H possesses VT symmetry 
(that is, that H commutes with the VT operator), and (ii) if (p is an eigenstate of H 
with eigenvalue E, then it is simultaneously an eigenstate of VT with eigenvalue A (it 
is this second assumption of unbroken VT symmetry that is crucial): 

H(l) = E(f) and VT(f) = X(f). (2) 

We first show that the eigenvalue A is a pure phase. We multiply VT(f) = Xcj) on the 
left by VT and use the fact that V and T commute and that = = 1 to conclude 
that (f) = and thus A = for some real a. Next, we introduce a convention used 
throughout this paper. Without loss of generality we replace the eigenstate by e~'"/^0 
so that its eigenvalue under the operator VT is unity: VTcp = 0. Next, we multiply 
the eigenvalue equation Hcj) = Ecf) on the left by VT and use [VT, if] = to obtain 
E(j) = E*(f). Hence, E = E* and the eigenvalue E is real. 

The crucial step in the argument above is the assumption that is simultaneously 
an eigenstate of H and VT. In quantum mechanics if a linear operator X commutes 
with the Hamiltonian H, then the eigenstates of H are also eigenstates of X. However, 
we emphasize that the operator VT is not linear (it is antilinear) and thus we must 
make the extra assumption that the VT symmetry of H is unbroken; that is, that is 
simultaneously an eigenstate of H and VT. This extra assumption is nontrivial because 
it is difficult to determine a priori whether the VT symmetry of a particular Hamiltonian 
H is broken or unbroken. For the Hamiltonian H in ^ the VT symmetry is unbroken 
when e > and it is broken when e < 0. The conventional Hermitian Hamiltonian for 
the quantum mechanical harmonic oscillator lies at the boundary of the unbroken and 
the broken regimes. Recently, Dorey et al. proved rigorously that the spectrum of H in 
is real and positive jl] in the region e > 0. Many other PT-symmetric Hamiltonians 
for which space-time reflection symmetry is not broken have been investigated, and the 
spectra of these Hamiltonians have also been shown to be real and positive (3]. 

It is important to know whether a given non-Hermitian "PT-symmetric Hamiltonian 
has a positive real spectrum, but the most urgent question is whether such a Hamiltonian 
defines a physical theory of quantum mechanics. By a physical theory we mean that 
there is a Hilbert space of state vectors and that this Hilbert space has an inner product 
with a positive norm. In quantum mechanics we interpret the norm of a state as a 
probability and this probability must be positive. Furthermore, we must show that the 
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time evolution of the theory is unitary. This means that as a state vector evolves in 
time the probability does not leak away. With these considerations in mind one would 
wonder whether a Hamiltonian such as if in (^J gives a consistent quantum theory. 
Indeed, early investigations of this Hamiltonian have shown that while the spectrum is 
entirely real and positive when e > 0, one inevitably encountered the severe problem of 
a Hilbert space endowed with an indefinite metric |BJ. 

However, there is a new symmetry that all PT-symmetric Hamiltonians having 
an unbroken PT-symmetry possess |7^. We denote the operator representing this 
symmetry by C because the properties of this operator resemble those of the charge 
conjugation operator in particle physics. This allows us to introduce an inner product 
structure associated with CVT conjugation for which the norms of quantum states are 
positive definite. Because of this we can say that VT symmetry is an alternative to 
conventional Hermiticity; it introduces the new concept of a dynamically determined 
inner product (one that is defined by the Hamiltonian itself). Consequently, we can 
extend the Hamiltonian and its eigenstates into the complex domain so that the 
associated eigenvalues are real and the underlying dynamics is unitary. This shows 
that PT-symmetric Hamiltonians are Hermitian in an extended (non-Dirac) sense. 

This paper is organized as follows. In Sec. 12 we give a general discussion of the C 
operator and in Sec. El we present a simple 2x2 matrix example of this operator. In 
Sec. m we show how to calculate C using perturbation theory for the cubic Hamiltonian 
H = + |/i^a;^ + iex^. In Sees. E] and IHl we calculate C for quantum mechanical 
Hamiltonians having two and three degrees of freedom. This is the principal new result 
in this paper. In Sees. [7|and|Hlwe consider possible physical applications and draw some 
conclusions. 



2. Construction of the C Operator 

We begin by summarizing the mathematical properties of the solution to the Sturm- 
Liouville differential equation eigenvalue problem 

- (p'^ix) + x'^{ixy(j)n{x) = En(t)n{x) (3) 

associated with the Hamiltonian H m This differential equation must be imposed 
on an infinite contour in the complex- a; plane. For large this contour lies in wedges 
placed symmetrically with respect to the imaginary-x axis j3] • The boundary conditions 
on the eigenf unctions are that exponentially rapidly as |x| — oo along the 

contour. For < e < 2, the contour may lie on the real axis. 

When e > 0, the Hamiltonian has an unbroken VT symmetry. Thus, we know that 
the eigenf unctions (l)n{,x) are simultaneously eigenstates of the VT operator: 

PT0„(a;) = A„0„(x). (4) 

As we argued above, A„ is a pure phase and, without loss of generality, for each n this 
phase can be absorbed into by a multiplicative rescaling so that the new eigenvalue 
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under VT is unity: 

PT0„(x) = = (5) 

It is not known rigorously yet, but there is strong evidence that when properly 
normalized the eigenfunctions (j)n{x) are complete. The coordinate-space statement of 
completeness reads 

= 6{x - y) (x, y real). (6) 

n 

This nontrivial result has been verified numerically to extremely high accuracy (twenty 
decimal places) jHUHl- There is a factor of (—1)" in the sum. This unusual factor does 
not appear in conventional quantum mechanics. The presence of this factor is explained 
in the following discussion of orthonormality [see (jH))] in which we encounter the problem 
associated with non-Hermitian PT-symmetric Hamiltonians. 

There seems to be a natural choice for the inner product of two functions f{x) and 
g{x): 

{f,g)^ J dx[VTf{xMx), (7) 

where VTf{x) = [f{—x)]* and the integration path is the appropriate contour in the 
complex-x plane. [We will see that (|7j) is not the correct choice for an inner product 
because it gives an indefinite metric. The correct inner product will be defined shortly, 
but studying this inner product is useful because it reveals the underlying mathematical 
structure of the theory.] The apparent advantage of this inner product is that the 
associated norm (/, /) is independent of the overall phase of f{x) and is conserved in 
time because H commutes with VT and the time-evolution operator is e~*^*. Phase 
independence is desired because in quantum mechanics the objective is to construct 
a space of rays to represent quantum mechanical states. With respect to this inner 
product, the eigenfunctions 0m(a;) and 0„(x) of if in (jT)) are orthogonal for n 7^ m 
because H is symmetric. However, when we set m = n we see by direct numerical 
calculation that the norm is evidently not positive: 

{<Pm, 0n) = i-irSmn. (8) 

This result is apparently true for all values of e in Q, and it has been verified numerically 
to extremely high precision. Because the norms of the eigenfunctions alternate in 
sign, the metric associated with the VT inner product (■, ■) is indefinite. This sign 
alternation appears to be a generic feature of this VT inner product. [Extensive 
numerical calculations verify that the formula in (jH)) holds for all e > 0.] We emphasize 
that while the sign of the norm of (pn is hard to verify analytically, the orthogonality 
of (pm and (pn is a trivial consequence of the symmetry of H. It is necessary to assume 
that H be symmetric in order to have this orthogonality. 

In spite of the nonpositivity of the inner product, it is instructive to proceed with 
the usual analysis that one would perform for any Sturm-Liouville problem of the form 
Hcpn = En(f>n- First, we use the inner product formula ^ to verify that ^ is the 
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representation of the unity operator. That is, we verify that J dyS{x — y)S{y — z) = 
6{x — z). 

Second, we show how to reconstruct the parity operator V in terms of the 
eigenstates. In coordinate space the parity operator is given by V{x,y) = 6{x + y), 
so from (jni) we get 

v{x,y) = J2i-^rM^)M-y)- (9) 

n 

By virtue of (jH)) the square of the parity operator is unity: = 1. 
Third, we reconstruct H in coordinate space: 

H{x, y) = 5^(-l)"i?„0„(a;)0„(y). (10) 



Using (jni) - (|H1) we can see that H satisfies H(j)n{x) = En(pn{x). 

Fourth, we construct the coordinate-space Green's function G{x,y): 

G{x, y) = ^(-l)"^0„(x)0„(y). (11) 
The Green's function is the functional inverse of H; that is, G satisfies 



dyH{x,y)G{y,z) 



+ x'^(ix) 



G{x,z) = 6{x - z). (12) 



dx"^ 

The time-independent Schrodinger equation cannot be solved analytically; it can 
only be solved numerically or perturbatively. However, the differential equation for 
G{x, z) in p2|) can be solved exactly and in closed form because it is a Bessel equation 
[H]. The technique is to consider the case < e < 2 so that we may treat x as real and 
then to decompose the x axis into two regions, x > z and x < z. We solve the differential 
equation in each region in terms of Bessel functions and patch the solutions together 
at a; = 2;. Then, using this coordinate-space representation of the Green's function, 
we construct an exact closed-form expression for the spectral zeta function (sum of the 
inverses of the energy eigenvalues). To do so we set z = x in G{x,z) and use (jH} to 
integrate over x. For all e > we obtain jH] 

r (ji;) r (4;) r (jfe) 



E- 

^ E 



^ , (1^) sin (^) 
cos (^) sin (1^) 



(4 + e)^r(i±f)r(|±f) 



(13) 



V4+e; 

This result has been verified to extremely high numerical accuracy [H]. 

All of these general Sturm-Liouville constructions are completely standard. But 
now we must address the crucial question of whether a PT-symmetric Hamiltonian 
defines a physically viable quantum mechanics or whether it merely provides an 
amusing Sturm-Liouville eigenvalue problem. The apparent difficulty with formulating 
a quantum theory is that the vector space of quantum states is spanned by energy 
eigenstates, of which half have VT norm +1 and half have "PT norm —1. Because the 
norm of the states carries a probabilistic interpretation in standard quantum theory, 
the existence of an indefinite metric in (jHJ seems to be a serious obstacle. The situation 
here in which half of the energy eigenstates have positive norm and half have negative 
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norm is analogous to the problem that Dirac encountered in formulating the spinor wave 
equation in relativistic quantum theory ^U]. Following Dirac's approach, we attack the 
problem of an indefinite norm by finding a physical interpretation for the negative norm 
states. We claim that in any theory having an unbroken VT symmetry there exists 
a symmetry of the Hamiltonian connected with the fact that there are equal numbers 
of positive-norm and negative-norm states. To describe this symmetry we construct a 
linear operator denoted by C and represented in position space as a sum over the energy 
eigenstates of the Hamiltonian [Zj: 

cix,y) = Y,Mx)My)- (14) 

n 

The properties of this new operator C closely resemble those of the charge 
conjugation operator in quantum field theory. For example, we can use (jHI) - (|H1) to 
verify that the square of C is unity (C^ = 1): 

dyC{x,y)C{y,z) = 5{x - z). (15) 

Thus, the eigenvalues of C are ±1. Also, C commutes with the Hamiltonian H. 
Therefore, since C is linear, the eigenstates of H have definite values of C. Specifically, 
if the energy eigenstates satisfy (jHl), then we have C0„ = (— 1)"0„ because 

C0n(a;) = / dyC{x,y)4)n{y) = ^(f)m{x) / dy (t)m{y)(pn{y)- (16) 

m 

We then use / dy (j)m{y)4>n{y) = {4>m, <Pn) according to our convention. We conclude that 
C is the operator observable that represents the measurement of the signature of the 
VT norm of a state. Note that since the C operator measures the VT norm of a state, 
we can think of the VT norm as the C "charge" of the state. 

The operators V and C are distinct square roots of the unity operator 5{x — y). 
That is, V^ = = 1, but V ^ C. Indeed, V is real, while C is complex. Note that the 
parity operator in coordinate space is explicitly real V{x,y) = 6{x + y); the operator 
C{x, y) is complex because it is a sum of products of complex functions, as we see in f|Tl|l . 
The complexity of the C operator can be seen explicitly in perturbative calculations of 
C{x,y) [TT]. We show how to perform these perturbative calculations in Sees. HI and El 
Furthermore, these two operators do not commute; in the position representation 

{CV){x,y) = ^(j)nix)(j)ni-y) but (PC)(x, y) = ^ 0„(-x)0„(?/), (17) 

n n 

which shows that CV = (VC)*. However, C does commute with VT. 

Finally, having obtained the operator C we define a new inner product structure 
having positive definite signature by 

{f\g)^ [ dx[CVTf{x)]g{x). (18) 

This inner product is phase independent and conserved in time like the VT inner product 
((7j). This is because the time evolution operator, just as in ordinary quantum mechanics, 
is e*^*. The fact that H commutes with the VT and the CVT operators implies that 
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both inner products, ((Tj) and (|18p. remain time independent as the states evolve in time. 
However, unhke ((Tj), the inner product (fTH|) is positive definite because C contributes 
— 1 when it acts on states with negative VT norm. In terms of the CVT conjugate, the 
completeness condition (jH)) reads 



Unlike the inner product of conventional quantum mechanics, the CVT inner product 

is dynamically determined; it depends implicitly on the Hamiltonian. 

The operator C does not exist as a distinct entity in conventional quantum 
mechanics. Indeed, if we allow the parameter e in (0) to tend to zero, the operator 
C in this limit becomes identical with V. Thus, in this limit the CVT operator becomes 
T, which is just complex conjugation. As a consequence, the inner product (fT8|l defined 
with respect to the CVT conjugation reduces to the complex conjugate inner product 
of conventional quantum mechanics when e — > 0. Similarly, in the e — limit, ()19p 
reduces to the usual statement of completeness J2n^n{x)(f)^{y) = 6{x — y). 

The CVT inner-product ()18|) is independent of the choice of integration contour C 
so long as C lies inside the asymptotic wedges associated with the boundary conditions 
for the Sturm-Liouville problem Path independence follows from Cauchy's theorem 
and the analyticity of the integrand. In ordinary quantum mechanics, where the positive- 
definite inner product has the form j dx f*{x)g{x), the integral must be taken along the 
real axis and the path of the integration cannot be deformed into the complex plane 
because the integrand is not analytic. [Note that if a function satisfies a linear ordinary 
differential equation, then the function is analytic wherever the coefficient functions of 
the differential equation are analytic. The Schrodinger equation is linear and its 
coefficients are analytic except for a branch cut at the origin; this branch cut can be 
taken to run up the imaginary axis. We can choose the integration contour for the inner 
product (jHI) so that it does not cross the positive imaginary axis. Path independence 
occurs because the integrand of the inner product (jH)) is a product of analytic functions.] 
The VT inner product (jZj) shares with (fTK|) the advantage of analyticity and path 
independence, but suffers from nonpositivity. We find it surprising that a positive- 
definite metric can be constructed using CVT conjugation without disturbing the path 
independence of the inner-product integral. 

Why are PT-symmetric theories unitary? Time evolution is determined by the 
operator e~*^* whether the theory is expressed in terms of a PT-symmetric Hamiltonian 
or just an ordinary Hermitian Hamiltonian. To establish the global unitarity of a theory 
we must show that as a state vector evolves, its norm does not change in time. If ipo{x) is 
a prescribed initial wave function belonging to the Hilbert space spanned by the energy 
eigenstates, then it evolves into the state i^tix) at time t according to 



With respect to the CVT inner product defined in (fTHjl . the norm of the vector ipt{x) 
does not change in time, (iptli^t) = (^olV^o); because the Hamiltonian H commutes with 




(19) 



n 



(20) 
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the CVT operator. 

Establishing unitarity at a local level is more subtle. Here, we must show that 
in coordinate space, there exists a local probability density that satisfies a continuity 
equation so that the probability does not leak away. This is a nontrivial consideration 
because the probability current flows about in the complex plane rather than along 
the real axis as in conventional Hermitian quantum mechanics. Preliminary numerical 
studies indeed indicate that the continuity equation is fulflUed [T^ . 

Just as states in the Schrodinger picture evolve in time according to the usual 
equation ((201), operators A also evolve according to the conventional Heisenberg-picture 
equation 

A{t) = e-'^U(O) e'^*. (21) 

Given this equation, it is clear how to deflne an observable in PT-symmetric quantum 
mechanics. The crucial property of an observable in any theory of quantum mechanics 
is that its expectation value in a state must be real. This will be true if 

_ jcvr ^ CVT A CVT, (22) 

where T signifles transpose. If this condition is fulfulled by a linear operator in a theory 
deflned by a PT-symmetric Hamiltonian, we say that the operator is "PCT symmetric" 
and is an observable. Note that this condition is the analog of the usual condition in 
conventional quantum mechanics that an observable be Hermitian in the usual Dirac 
sense A^ = A: 

A^ = A*. (23) 

The condition for an operator to remain an operator as time evolves is simply that H 
be symmetric: = H. This symmetry condition has been implicitly assumed in all 
of the "PT-symmetric models discussed in the literature. Recall that the symmetry of 
H ensures that eigenstates of H corresponding to different energies will be orthogonal. 
Note that there are two time-independent observables in the theory; namely, H and C. 

3. Illustrative Example: A 2 x 2 Matrix Hamiltonian 

We illustrate the above results concerning PT-symmetric quantum mechanics using the 
finite-dimensional 2x2 symmetric matrix Hamiltonian 

where the three parameters r, s, and 6 are real. This Hamiltonian is not Hermitian in 
the usual Dirac sense, but it is VT symmetric, where the parity operator is ^iHj 



: (25) 




and T performs complex conjugation. (Note that T does not perform Hermitian 
conjugation, or else it would not leave the commutation relation [x,p] = i invariant.) 



The C Operator in VT -Symmetric Theories 



11 



There are two parametric regions for this Hamiltonian. When < sin^ 6, 
the energy eigenvalues form a complex conjugate pair. This is the region of broken 
VT symmetry. On the other hand, if > sin^ 9, then the eigenvalues e± = 
r cos — sin 9 are real. This is the region of unbroken VT symmetry. In 

the unbroken region the simultaneous eigenstates of the operators H and VT are 

if e^"/2 \ , , , if e-*"/2 

= /o P-W2 1^-) = /o 

V 2 COS a \ e ' / v 2 cos a 

where we set sin a = {r / s) sm9. It is easily verified that {s±,e±) = ±1 and that 
{e±,e^) = 0, recalling that {u,v) = (VTu) ■ v. Therefore, with respect to the VT 
inner product, the resulting vector space spanned by energy eigenstates has a metric of 
signature (+, — ). The condition > sin^ 9 ensures that VT symmetry is not broken. 
If this condition is violated, the states (pUj) are no longer eigenstates of VT because a 
becomes imaginary. (When VT symmetry is broken, we find that the VT norm of the 
energy eigenstate vanishes.) 

Next, we construct the operator C: 




C = . (27) 

COS a 

Note that C is distinct from H and V and has the key property that C\6±) = ±\e±). The 
operator C commutes with H and satisfies = 1. The eigenvalues of C are precisely 
the signs of the VT norms of the corresponding eigenstates. Using C we construct the 
new inner product structure {u\v) = {CVTu) ■ v. This inner product is positive definite 
because {e±\e±) = 1. Thus, the two-dimensional Hilbert space spanned by \e±), with 
inner product (-I-), has a Hermitian structure with signature (+, +). 

Let us demonstrate explicitly that the CVT norm of any vector is positive. For 
the arbitrary vector ip = (°), where a and b are any complex numbers, we see 
that T^ = (ll), that VTtP = (JI), and that CVTtP = ^ Hl-Zt^a)- Thus, 
{ip\ip) = (CVTip) ■ i/j = ■^^[a*a + b*b + i{b*b - a* a) sin a]. Now let a = x + iy and 
b = u + iv, where x, y, u, and v are real. Then 

{ip\ip) = (x^ + f ^ + 2xv sin a + y'^ + — 2yu sin a) / cos(q;), (28) 

which is explicitly positive and vanishes only i{x = y = u = v = 0. 

Since {u\ denotes the CVT conjugate of \u), the completeness condition reads 

\e^){e^\ + \e-){s-\=(^l J (29) 

Furthermore, using the CVT conjugate {e±\, we have C as C = |£:+)(£+| — \e^){e^\, 
instead of the representation in (fT^ . which uses the VT conjugate. 

For the two-state system discussed here, if —>■ 0, then the Hamiltonian ()24p 
becomes Hermitian. However, in this limit C reduces to the parity operator V. As a 
consequence, the requirement of CVT invariance reduces to the standard condition of 
Hermiticity for a symmetric matrix; namely, that H = H* . This is why the hidden 
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symmetry C was not noticed previously. The operator C emerges only when we extend 
a real symmetric Hamiltonian into the complex domain. 

4. Perturbative Calculation of the C Operator for an ix^ Theory 

The C operator can be calculated in some infinite-dimensional quantum mechanical 
models. For an + ix^ potential C can be obtained from the summation in (fT^ 
using perturbative methods and for an — potential C can be calculated using 
nonperturbative methods fT]. In this paper we focus on perturbative methods for 
calculating C. 

Let us consider the "PT-symmetric Hamiltonian for a harmonic oscillator perturbed 
by an imaginary cubic potential: 

H = \p^ + ix^ + iex^ (e real). (30) 

Following the procedure in Ref. jTT], we note that the energy eigenstates are 
solutions of the Schrodinger equation 

H(t)n{x) = (^-l^ + + ^^^^ 0n(a;) = E^M^)- (31) 

The eigenstates and the corresponding eigenvalues may be expressed as series in powers 
of e by perturbing around the known energy eigenstates and eigenvalues of the harmonic 
oscillator. To second order in perturbation theory the eigenstates have the form 

M^) = e--^""" K(x) + eP„(x) + e'Q^x)] (32) 

with energies given by 

E^ = n + ^ + eA^ + e^B^. (33) 

Here, Pn{x) and Qn{x) are polynomials in x of degree n + 3 and n + 6, respectively; 
a„ is a normalization constant to be determined. We include a factor of because the 
unperturbed wavefunctions have the form 



.1^.2 



= _^ e-t--ff,(x), (34) 

where Hn{x) are Hermite polynomials. This ensures that the unperturbed wavefunctions 
are eigenstates of the VT operator with unit eigenvalue: 

pr0f(x) = 0f(x). 

The wave functions are PT-normalized according to 

dx [<p^:\x)]' = {-ir. 
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4-1. First- Order Calculation of the Energy Eigenstates and Eigenvalues 
Note that to order e°, (iBlll becomes 
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-1^ (i^ d 
^ dx"^ dx 

and to order e\ we have 

1 rf2 



n 



Hn{x) 



(35) 



_l_ — ri — - 

^dx^^ 2 



e 2 



ixne 2^' if„(x) 



or 



1 d^ 



d 



n 



PJx) 



Since the Hermite polynomials form a complete set, we may rewrite any polynomial 
as a linear superposition of Hermite polynomials. Rewriting Pn{x) in this manner yields 



= {An - ix^)Hn{x), 



which with the help of ()35|) simplifies to 

"^Pkik - n)Hk{x) = {An - ix^)Hn{x). 



(36) 



Also, we have 



x^Hnix) 



13 3 

-Hr,4-?.{x) + -(n + l)Hn+i{x) + -n'^Hn-i{x) 



+ ri(n-l)(n-2)/J„_3(x). (37) 

The coefficient of Hn{x) on the left side of equation (jHUj) is zero, and the expression for 
x^Hn{x) on the right side does not contain any terms in Hn{x). Hence, we conclude 
that An = for all n to first order in perturbation theory. Thus, the perturbed energy 
equals the unperturbed energy. 
Rewriting equation as 



^Pkik - n)Hk{x) 



^Hn+z{x) + -^{n + l)Hn+i{x) 



+ ^n^Hn^i{x) + n{n - l){n - 2)Hn^-i{x)\ (3^ 



and comparing coefficients reveals that 



1 3 

iPn{x) = —Hn+six) + -{n + l)Hn+lix) 

3 1 

- -n^Hn~i{x) - -n{n -l){n- 2)iJ„_3(a;) 



(39) 
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4-2. Second- Order Calculation of the Energy Eigenstates and Eigenvalues 
At order e^, the eigenproblem becomes 
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d^ 



1_ — n — - 

^dx^^ 2 



or 



1 d^ 



d 



^ dx"^ ~'~ dx 



Qn{x) = BnHn{x) - ix^ Pn{x) . 

On posing Qn{x) = J2k1kHk{x) this reduces to 

^qk{k - n)Hk{x) = BnHn{x) - ix^Pn{x). 



(40) 



Combining equations (jHTj) and (jniJ|) we obtain 



+ i (7^^ + 33n + 27) if„+2(x) + ^ (SOn^ + 30n + ll) 
lb 8 

- - 1) {7n^ - 19n + l) Hn-2{x) 
1 

- -n(n - l)(n - 2)(n - 3)(4n - ?>)Hn^i{x) 

- \{n - l)(n - 2){n - 3)(n - 4)(n - 5)//„„6(x). 



(41) 



Substituting this result back into (pOj) and comparing coefficients, we find that 

Qn{x) = 



-^—Hn+eix) —(An + 7)Hn+4{x) 

1152 ^ ^128 



^ (7n2 + 33n + 27)ff„+2(x) 



n(n 



1) (7^2 - 19n + 1) /f„_2( 



a; 



- -n{n - l)(n - 2)(n - 3)(4n - 3)Hn-4{x) 
8 

1 



l; 



-n{n -l){n- 2){n - 3)(n - 4)(n - 5)i/„_6(a;), 



(42) 



and 



E„ 



^+1 + ^(SOn^ + 30n + 11)6^ + O(e^). 
2 8 



Having found expressions for the eigenfunctions of the Hamiltonian, we must verify 
that they are PT-normahzed to this order in perturbation theory: J^^dx [(f)n{x)]'^ = 
(—1)" + C(e^). This determines the value of a„ in 

^2 



2"n! 



1. 



2"n!0F V 7-00 ^ ^ 

Using ()39|) as well as the orthogonality and normalization conditions for the Hermite 
functions, we obtain 1 = [l — {^n^ + + ^-n + ||) e^] , and hence 

a„ = 1 + Y^(2ri + l)(82n2 + 82n + ^7)e^ + C(e^). 
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We must also verify that the energy eigenstates are simultaneously eigenstates of 
VT. Note that Hn{x) and Qn{x) are even in x for even n, and odd in x for odd n; 
Pn{x) has the opposite parity, but it contains an additional factor of i. Hence, all three 
polynomials are PT-symmetric for even n, and VT anti-symmetric for odd n. 



VTH^{x) 



(43) 



Hn{x) if n is even, 
—Hn{x) if n is odd. 

The same holds for the prefactor z", and thus (f)n{x) is indeed PT-symmetric for all n 
VT(t)n{x) = (f)nix). 



4-3. Calculation of the C Operator 

We can now construct the C operator for the ix^ theory to order 0{e): 

oo 

C{x,y) = 



n=0 



1 °° ( — 1)'^ 



n=0 
oo 



2"n! 



1 e-^(^^+^^)5^-^ff„(x)/7„(-y) 

V ^— n 



+ e 



n=0 

oo 



e-^(^^+^^) J2 :^Pn{x)H^{-y) + {x^y) 



n=0 



C(e2). (44) 



To proceed, we must make use of the completeness relation for the Hermite functions. 
To that end, we first need to express Pn{x) solely in terms of Hn{x) and derivatives 
thereof. A comparison of equations (j37j) and (j39|) shows that 

1 1 

iPn{x) = —Hn+3{x) + Hn+i{x)--x^Hnix)+2xnHn{x)-3n^Hn-i{x).iA5) 

We now use Hn+i{x) = (2x — ^) Hn{x) to obtain 



Also, from 2nHn{x) = 



d di , o 1 d) , '~> 1 

— + 6x— + 6(1 - 2x^)— + A{2x^ - 3x) 
dx-^ dx^ dx 



HJx). 



dx 



Hn{x), we have 



-3n^Hn-i{x) = --nH'^{x) 



3 d^ 3 d^ 



3 d 



-X 



Adx^ 2 dx"^ 2dx 



Substituting all this into (pSj) gives 
iPJx) ■- 



2 d^ 2x^^+(x2 2)^+(^x^ + x 

3 dx^ dx"^ dx \ 3 



Hn{x). 



HJx) 
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1^2 I 2 



-TT + 32;-^ + 3(1 + x^)-^ + {x^ + 3x) 
dx-^ dx^ dx 



- 2x ( + 2x-^ + (1 + x"^) 
ax^ dx 



+ {x' 



X + 



d ' 

dx 



+ \ X + -a;' 



TT 



Finally, we obtain an expression of the form 

xPri \ X 



62"^' 5. ( e ^''^HJx] 



X 



where the differential operator 6^ is given by 
'2 d^ _ 2d 

3 dx^ dx 

It is now easy to complete our calculation of the C operator. We simply 
our new expression for Pn{x) into equation (jUj): 



C{x,y) 



4- 4- Verification of C 



Ad' d 
3 dx-^ dx 



5{x + y). 



[x) j . (46) 

(47) 

(48) 
substitute 

(49) 



We can perform several checks to ascertain the correctness of the C operator in (|49|) . 
For example, we can verify that the eigenstates of H are also eigenstates of C: 

(50) 



dyC{x,y)(t>n{y) = {-lT(t>n{x). 
Now, equation (fSnjl is clearly satisfied to zeroth order in e: 

dy5{x + y)e--^y"Hr,{y) = {-1)^6" H^x) 



because if„(— x) = (— l)"if„(x). 
To first order in e, ()50p reads 



-l)"e-^-'R,(x 



dyS{x + y)e P„(?/) 

) 



Noting that P„(— x) = — (— l)'^P„(x) we can write 



2i(-l)"e-5^'P„(x) 



or 



iPJx) 



1„2 
62^ 



4 /""^ 1 

-x^ / dy6{x + y)e-^y'Hn{y) 

1 

+ 2xpx / dy 5{x + y)ye~2y Hn{y) 
- xpx (xe~^^^ Hn{x) 



(51) 

(52) 
(53) 



^x^ (^e 2'^'if„(x 



which agrees with the previous result in (j47p . 
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4-5. The C Operator as an Exponential 

Extending the calculation of the C operator to second order in perturbation theory 
presents no new conceptual difficulties over and above the ones encountered at first 
order in perturbation theory. We simply cite the result given in jTT]: 

C{x,y) = 1 — e ( -p"^ — 2xyp 

+ (^^p' - ^xyp' + i2xY - 12)/J J 5(x + y) + ^(e^), (54) 

where p = —i-^- The structure of this formula suggests that it might be rewritten as 
an exponential: 



C{x, y) = exp 



^4 , \ 
\-p— 2xpx j 



6{x + y). (55) 



Observe that the C operator reduces to V in the limit where the parameter e tends to 
zero in ()49p . Note that the expression for the parity operator V is V = e*™^"^, where 
and a represent the standard quantum mechanical harmonic oscillator raising and 
lowering operators, respectively. The combination a^a represents the number operator. 
It is interesting that in exponential form the C operator is a series in odd powers of e. 
(Consult Ref. [TT] for the contribution to C) 

5. Perturbative Calculation of the C Operator for an ix^y Theory 

Let us now apply the techniques of the previous section to a quantum mechanical theory 
having two degrees of freedom. Again, the PT-symmetric Hamiltonian consists of a 
standard harmonic oscillator interacting with a complex potential 

^ dx"^ ^ dy"^ 

This complex Henon-Heiles theory was studied in Ref. 

The chain of reasoning culminating in an explicit expression for the C operator 
is much the same as that in the previous section, allowing us to concentrate on the 
essentials. In order to solve the Schrodinger equation 

H(j)mn{x, y) = Emn(f)mn{x, y) , (57) 

we again resort to perturbative methods. To ffist order in e we have 

(J^mnix, y) ~ e-^(^-'+^') [H^{x)H^{y) + eP^„(x, y)] (58) 

and 

Emn = m + n+l + eAmn- (59) 



H = -k^-kTrT + kx' + h +i^x'y. (56) 
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1 92 1 92 d d 

+ X— — h y— m — n 



2 dx"^ 2 dy"^ dx dy 
Rewriting Pmn{x,y) as a power series in Hermite polynomials 

Pmnix,y) = ^PrsHr{x)Hs{y), 



Pmn{x, y) = {Amn " ix^ y) H m{x) H n{y) ■ 



then yields 



^Prs [{r - m) + (s - n)] Hr{x)Hs{y) = (A„„ - ix^y)Hm{x)Hn{y). (60) 



We then have 

1 1 

x'^yHm{x)Hn{y) = -Hm+2{x)Hn+l{y) + -nHm+2{x)Hn-l{y) 

+ ^ + H^{x)Hn+i{y) + + nHm{x)Hn-i{y) 

+ ^m{m - l)Hm~2{x)Hn+i{y) + m(m - l)nHrri-2{x)Hn-i{y). (61) 

The right side of this equation does not contain any terms in Hm{x)Hn{y), allowing us 
to deduce that A^n = for all m and n. Thus, the energy does not change to first order 
in perturbation theory. 

A comparison of the coefficients in (jHUI) reveals that 



iPmn{x,y) = Hm+2{x) 



■ 1 1 

—Hn+i{y) + -nHn-i{y) 



+ I m + - ) Hm{x) 



-Hn+i{y) - nHn-i{y) 



- m{m - l)Hm-2{x 
We can rewrite this equation in the form 

iP^nix.y) = e^(^'+^')5., [e-'^^^'+y'^H^{x)HM 

where the differential operator 5xy is given by 

2 93 1,9291 



]^Hn+i{y) + ^nHn^i{y) 



1 2 d 
-X — xy— y. 



(62) 



(63) 



(64) 



Sdx'^dy 3"" dy 3"^^ dx 3 
5.2. Calculation of the C Operator 

Having established the form of equation (jUHj) . it is now straightforward to calculate the 
C operator: 

00 00 



C{x, x'; y, y') = ^ y)(t'mn{x', y') 



m=0 71=0 
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[1 - ie {5^y + S^>y>)] 5{x + x')S{y + y') 



4 ^3 
1 - ^e - 



d A d 



Sdx^dy ~'~ 3 dy 3 ^ dx 
where the parity operator is 6{x + x')6{y + y'). 



Six + x')Siy + y'), 



(65) 



6. Perturbative Calculation of the C Operator for an ixyz Theory 

As a third example we consider a quantum mechanical theory with three degrees of 
freedom. Once again, the PT-symmetric Hamiltonian consists of a standard harmonic 
oscillator part interacting with a complex potential 

id' 



H 



^ dx'^ ^ dy'^ 



We wish to solve the Schrodinger equation 
H4>kim{x, y, z) = Ekim4>kim{x, y, z), 
whose eigenfunctions are given by 

H,{x)Hi{y)Hm{z) 
+ (^Pkimix, y, z) + e'Qkim{x, y, z) 



1 2 

kimix^y, z) ~ e"2'' 



and whose energies have the form 

3 
2 

where we have set = + j/^ + z"^ and n = k + I + m. 



(66) 
(67) 



(68) 



6.1. Energy Eigenstates 

To first order in e (|67p becomes 

1 92 1 92 , 92 d d d 



+ X— + y— + z- n Pkimix, y, z 



^ dx'^ ^ dy'^ ^ dz'^ ' dx " dy ' dz 
= {Akim - ixyz)Hk{x)Hi{y)Hm{z). 

Rewriting Pkim{x, y, z) as a sum of Hermite polynomials, 

Pkim{x,y,z) = ^prstHr{x)Hs{y)Ht{z), 



(69) 



r,s,t 



we obtain the equation 

"^Prstir + s + t- n)Hr{x)Hs{y)Ht{z) 



r,s,t 



{Akim - ixyz)Hk{x)Hi{y)Hm{z). 



(70) 
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1 

+ 4 



kHk-i{x)Hi+i{y)Hm+i{z) + lHk+i{x)Hi^i{y)Hm+i{z) 



+ mHk+i{x)Hi+i{y)Hm^i{z) 



1 

+ 2 



+ kmHk-i{x)Hi+i{y)Hm-i{z) + lmHk+i{x)Hi^i{y)Hjn-i{z) 

+ klmHk-i{x)Hi^i{y)H,n-i{z). (71) 

The right side of this equation, being devoid of terms in Hk{x) Hi{y) Hm{z) , confirms that 
the energy is unaltered to this order in perturbation theory. In other words, Akim = 
for all k,l,m. Comparing coefficients, we then find that 

iPkim{x,y,z) = ^Hk+i{x)Hi+i{y)Hm+i{z) 



1 

+ 4 



kHk^i{x)Hi+i{y)Hm+i{z) + lHk+i{x)Hi^i{y)Hm+i{z) 

+ mHk+i{x)Hi+i{y)Hm^i{z) klHk^i{x)Hi_i{y)Hm+i{ 
+ kmHk-iix)Hi+i{y)Hm-iiz) + lmHk+i{x)Hi_i{y)H^_i{z) 



- -klmHk^i{x)Hi^i{y)R^^i{z). 



(72) 



Having established the form of Pkim{x,y,z), we now must ensure that the 
wavefunctions (j68|) are correctly PT-normalized to order e in the sense that 




dx dy dz [(f)kimix,y, z)y 



(73) 



Note that Pkim{xyz) does not contain a term in Hk{x)Hi{y)Hm{z). Hence, from the 
orthogonality and normalization conditions for the Hermite functions, it follows that 
the correctly normalized wavefunction must take the form 



')kim[x,y, Z) 



^'K^m^kmm. 



■.e [Hk{x)Hi{y)Hra{z) + ePkim{x,y,z)] 



6.2. The C Operator 

Once again, we must express Pkim{x,y, z) solely in terms of Hk{x)Hi{y)Hm{z) and 
derivatives thereof before we can apply the standard completeness relation for the 
Hermite functions. The result is 



iPkim{x,y,z) 



Q3 



3 dxdydz 3 \ dydz ^ ^ dxdz ^ dxdy ) 



d 



d 



d 



+ - I xy— — h xz— — h yz— 



dz 



dy 



dx 



-xyz 



Hk{x)Hi{y)Hm{z). (74) 
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An equivalent form of the polynomial, but one that is more amenable to the task 
in hand, may be obtained from 

d " 



nHn-i{x) = ^6 2^^" 



Hn+i{x) = 

ax \ \ 

We can thus write iPkimix, y, z) in the compact form 

iPkim{x,y,z) = e^'''6^y, e~^''^ Hk{x)Hi{y)H^{z) 



(75) 



(76) 



where the differential operator 5xyz is given by 



2 



d 



d 



d 



^^^^ 3 dxdydz 3 dz ^ dy ~'~ dx 

Given the symmetric nature of the operator 6xyz it is now particularly easy to derive 
the form of the C operator: 

OO CO oo 

C{x, x'; y, y'; z,z') = <Pkim{,x, y, z)(f)kimix', y', z!) 

k=0 1=0 m=0 

= [l-ie (S^y^ + 5x>y'^>)]S{x + x')5{y + y')5{z + z') 
= (1 - 2ie6,y,) 6{x + x')6{y + y')6{z + z') 

A 2/(9 d d 



1 -ie 



3 dxdydz ?>\^dz^ dy^^dx 



5{x + x')5{y + y')5{z + z').{n) 



6.3. Verification of C 



It is important to verify that the eigenstates of the Hamiltonian are also eigenstates of 
the C operator: 



(-l)"0Hm(a;, y,z) = JJJ dx' dy' dz'C{x, x'; y, y'; z, z')4>kim{x', y', z!). (78) 
To demonstrate this to first order in e, we have 

(-l)"e"5'''Pfc/m(a;, y-,z)= [ [ [ dx' dy' dz' 5x,x'5y,y'5z,z'e~^'''^ Pkim{x' , y', z') 

(79) 




-2iSxyz dx' dy' dz' 6x,x'Sy,y'6z^^'e ^""'^ Hk{x')Hi{y')Hm{z'), 



where we have used the abbreviated notation S^^x' ioT 6{x + x'). 

Observing that Pkimi-x, -y, -z) = Pkim{x, y, z), we can write 



2,i6xyz 



-ire—' Hk{x)Hi{y)H^{z) 



2{-ire-t' PkiM.y.z) 
whence we recover equation (f7(T|) : 

iPkim{x, y, z) = e5'''5(x, y, z) e'^''^ Hk{x)Hi{y)H^{z 
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6.4- Second Order Perturbation Theory - Difficulties with Degeneracy 

At order a tough problem surfaces, namely that of degeneracy. To second order the 
eigenproblem (|H7jl becomes 

J2 "Irstir + s + t- n)Hr{x)H,{y)Ht{z) 

r,s,t 

= BnHk{x)Hi{y)Hm{z) - ixyzPkim{x, y, z), (80) 

where we posed 

Qkim{x, z) = ^ qrstHr{x)Hs{y)Ht{z). (81) 

In order to find the coefficients grst, we need to express ixyzPkim{x,y, z) in terms 
of Hermite polynomials. We use the formula 

xHn+l{x) = lHn+2{x) + (u + l)if„(x) 

applied to (ff^ . The result is a highly symmetric formula that is too long to give here. 

When we now examine (jHUj) in the light of the formula for ixyzPkim{x,y, z), it 
soon becomes apparent that we run into an unforseen problem: the left-hand side 
of (jHUj) is zero whenever r + s + t = n, but at the same time we have terms like 
k{k — 1) Hk-2{x) Hi{y) Hm+2{z) and permutations thereof on the right-hand side which 
are clearly not zero. The underlying cause of this mismatch lies in the symmetric 
nature of the Hamiltonian and the associated degeneracy of its eigenvalues. In fact, the 
unperturbed eigenvalues En^ = n + 3/2 are {n + l){n + 2)/2-fold degenerate, where n 
denotes the energy level. While this degeneracy persists to first order in perturbation 
theory, it is partially lifted at second order. As a result of this degeneracy, one needs to 
take into consideration a mixing of states corresponding to the same energy. 

We briefly illustrate the technique here by examining the energy levels n = 2, 4, 
and 6, which are 6, 15, and 28 fold degenerate, respectively. A little reflection reveals, 
however, that not all of these states figure in the mixing. For the n = 4 level, for 
example, we find that we need only include 6 of the 15 states in the mixing. In the 
study of our three examples we shall need to have recourse to some special cases of the 
lengthy equation for ixyzPkim{x, y, z). 

6. 5. The n = 2 Energy Level 

We can repair the inconsistency encountered in ()80|) by replacing it here with 
J2 Qrst [{r + s + t- 2)] Hr{x)Hs{y)Ht{z) 

r,s,t 

= B2[a,H2ix)Hoiy)Hoiz) + a2H^{x)H2{y)H^{z) + aM^)Ho{y)H2{z)] 

- ixyz[aiP2m{x, y, z) + 02^020(2;, V-, z) + a3Poo2(a;, z)]. (82) 

We need to choose the mixing coefficients ai, a2, and 03 such that the problem terms 
H2{x)Ho{y)Ho{z), HQ{x)H2{y)Ho{z), and Ho{x)Hq{ii)H2{z) disappear. This amounts 
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(83) 



-1/4 B2 - 3/8 
y -1/4 -1/4 B2 

Cramar's rule states that for a nontrivial solution to exist we must require that the 
determinant of the given matrix be zero: 

{B2 - 1/8)2(52 - 7/8) = 0. 

We see that the effect of the second order contribution to the unperturbed energy level 
(of value 3.5 and 3- fold degenerate with respect to the states under consideration) is to 
split it into two levels, both of which are raised and one of which is doubly degenerate. 
To the doubly degenerate energy there corresponds the following condition on the mixing 
coefficients: 

ai + ^2 + as = 0, (e.g. ai = 1, 0^2 = 0:3 = -1/2). 

This reflects an anti-symmetric mixing of states. For the nondegenerate energy, the 
condition reads: 



tti = ^2 = as, 



(e.g. ai = a2 = as 



indicating a symmetric mixing of states. 

Let us consider the case of the symmetric mixing of states for illustrative purposes. 
Equation (|H^ becomes 

Qrst [{r + s + t- 2)] Hr{x)Hs{y)Ht{z) 



r,s,t 



= l[H2{x)Ho{y)Ho{z) + Ho{x)H2{y)Ho{z) + Ho{x)Ho{y)H2{z)] 

o 

- ixyz[P2oo{x, y, z) + Po2o(a;, y, z) + ^002(3;, y, z)]. 
From ()72|) we know that 

^[^200(2;, y, z) + Po2o(a;, y, z) + Poo2(a;, y, z)] 
= h^H^{x)HMH^{z) + H^{x)HMH^{z) + H^{x)HMH-i{z)] 

+ ^-H,{x)H,{y)H,{z), 
and in addition, we have 

ixyz [P2m{x, y, z) + Po2o(a;, z) + Poo2(a;, y, z)\ 
\^^{x)U2{y)U2{z) + H^{x)H2{y)H^{z) + H2{x)H2{y)H^{z 

H^{x)Ho{y)Ho{z) + Ho{x)H4{y)Ho{z) + Ho{x)Ho{y)H,{z 

H2{x)H2{y)Ho{z) + H2{x)Ho{y)H2{z) + Ho{x)H2{y)H2{z) 

H2{x)Ho{y)Ho{z) + Ho{x)H2{y)Ho{z) + Ho{x)Ho{y)H2{z) 



15) 



192 
1 

+ 48 

1 

+ 2L 
7 

+ 77 
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+ - H^{x)H2{y)Ho{z) + H^{x)Ho{y)H2{z) 
9o L 

+ H2{x)H^{y)Ho{z) + H2{x)Ho{y)H^{z) 
+ Ho{x)H^{y)H2{z) + Ho{x)H2{y)H,{z) 

+ ^H,{x)H2{y)H,{z) + ^H,{x)H,{y)Ho{z). 



(86) 



Having done this analysis, we now observe that the problem terms H2{x)Ho{y)Ho{z), 
HQ{x)H2{y)HQ{z), and Hq{x)Hq{ii)H2{z) do indeed disappear. 

Finally, we can determine the coefficients g^st in ()82|1 . We find that 

Q2oo(a;, z) + Qo2o(a;, 2;) + (5002(2;, y, z) 

1 



1152 



[H^{x)H2{y)H2{z) + H2{x)H,{y)H2{z) + H2{x)H2{y)Hi{z)] 



- -[H,{x)Hoiy)Ho{z) + Ho{x)H^{y)H^{z) + Ho{x)Ho{y)H^{z)] 
96 

- ^[H2{x)H2{y)Ho{z) + H2{x)Ho{y)H2{z) + Ho{x)H2{y)H2{z)] 



1 
384 



H,{x)H2{y)Hoiz) + H^{x)H^{y)H2{z) 
+ H2{x)H^{y)Ho{z) + H2{x)Ho{y)H,{z) 
+ Ho{x)H,{y)H2{z) + H^{x)H2{y)H,{z) 

- ^H2{x)H2{y)H2{z) + ^ifo(x)i/o(l/)i^o(^). 
(J. (J. r/ie n = 4 Energy Level 

At this energy level we assume a structure of the form 

qrst [(r + s + t - 4)] Hr{x)Hs{y)Ht{z) 

r.s.t 



(87) 



= 54 [«ii/4(a;)i/o(?/)i^o(2:) + a2Ho{x)H^{y)Ho{z) + aMx)Ho{y)H^{ 
+ (3Mx)H2{y)H2{z) + l32H2{x)H,{y)H2{z) + P,H2{x)H2{y)Ho{z 



— ixyz 



aiPiOQ{x, y, z) + 02^040 (a;, 2:) + a3Poo4(a;, y, z) 



(88) 



/?iPol22(a;, y, z) + /?2-P202(a;, y, 2;) + /?3-P22o(a;, ^) 

Now, we need to choose the mixing coefficients ai, 0:2, as, /32, and /^s such that 
the six problem terms disappear. Equivalently, we need to solve the system of linear 
equations M.v = 0, where the coefficient matrix M of the system is given by 

Ba-]^ 

-i 



M 



^4-24 u u -\ 

i?4-i -\ -\ 

54-^ -T -T 





_ 17 _1 _1 

4 24 4 4 

3 _3 R _ 11 _1 _i 

2 2^8 4 4 

I -I -\ B.-^ -\ 



(89) 
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and V denotes the column vector (ai, a2, a^, Pi, P2, Ps)- 

For a nontrivial solution we require that the determinant of the given matrix be 
zero, so 

(192S| - 496^4 - 33)(192S| - 352^4 + 81)^ = 0. 

The unperturbed energy level (of value 5.5 and 6-fold degenerate) is split into four 
distinct levels (3 raised, 1 lowered), of which two are doubly degenerate. 

Associated with the doubly degenerate values of the energy are the following 
conditions on the mixing coefficients: 

ai + q;2 + q;3 = and /3i + + /^g = 0, with /3i = ^(-5 T V2Ai)ai. 

b 

Hence, this case corresponds to an anti-symmetric mixing of states. 
For the nondegenerate energies, one has 

Oil — a2 — Ois and (3i — P2 — P3, with Pi — -{7 ± \/265)q;i, 

which yields a symmetric mixing of states. We observe that the unperturbed energy 
level (of value 7.5 and 10-fold degenerate) is split it into seven levels (6 raised, 1 lowered), 
of which three are doubly degenerate. 

To i^e = I there correspond the conditions 

ai = 0:2 = as = 0, /3i = /33 = Pb, Pi^ Pe, 

Pl + P2+P3 + Pi + P5+P6 = 0, -f = 0. 

We find that an odd permutations of the indices of the states introduces a relative minus 
sign. 

The symmetry of the states associated with the nondegenerate energies {Bq — 
5.473, 2.343, 0.391) is characterized by 

0/^ = 0/2 = as, Pi^ P2^ Ps^ Pi^ P^^ Pe, 

with the relations 

Pi = 8.863ai and 7 = 3.701/3i = 32.804ai, 

pi = 2.603ai and 7 = -12.888/5i = -33.553ai, 

/5i = -1.300ai and 7 = -3.396/5i = 4.415ai, 

respectively. Finally, for the case n = 6 the states corresponding to the doubly 
degenerate energies {Bq = 4.003, 1.981, —0.193) mix according to 

ai + a2 + a3 = 0, Pi + P2 + p3 + p4 + p5 + Pa = 0, 7 = 0. 

On the basis of the three cases considered, one can see that there is a direct 
correlation between the degeneracy of an energy level and the mixing symmetry of its 
associated state. An anti-symmetric mixing of states corresponds to a doubly degenerate 
eigenvalue, while a symmetric mixing of states or a state of mixed symmetry corresponds 
to a nondegenerate energy. Moreover, the number of (not necessarily distinct) energies 
equals the number of states we are mixing (in our cases: 3, 6, and 10). 
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In conclusion, it is necessary to take great care to deal with the difficulties presented 
by degeneracies. Through the examination of three examples we have found that the 
states corresponding to a degenerate energy mix according to certain symmetry criteria. 
These clearly have to be respected when one is attempting to calculate the C operator. 
Obviously, the problem of degenerate states makes it very difficult to calculate the C 
operator in systems having more that one degree of freedom. The problems associated 
with degeneracy can, in fact, be overcome and the techniques for doing so are described 
in a paper under preparation by Bender, Brody, and Jones jTH]; in this paper it is shown 
that it is even possible to find C for systems having an infinite number of degrees of 
freedom (quantum field theory). 

7. Applications and Possible Observable Consequences 

We do not know if non-Hermitian, PT-symmetric Hamiltonians can be used to describe 
experimentally observable phenomena. However, non-Hermitian Hamiltonians have 
already been used to describe interacting systems. For example, Wu showed that 
the ground state of a Bose system of hard spheres is described by a non-Hermitian 
Hamiltonian ^Hj- Wu found that the ground-state energy of this system is real and 
conjectured that all energy levels were real. Hollowood showed that even though the 
Hamiltonian of a complex Toda lattice is non-Hermitian, its energy levels are real |17j . 
Non-Hermitian Hamiltonians of the form H = + ix^ also arise in Reggeon field theory 
models that exhibit real positive spectra ^H] • In each of these cases the fact that a non- 
Hermitian Hamiltonian had a real spectrum appeared mysterious at the time, but now 
the explanation is simple: In each case the non-Hermitian Hamiltonian is PT-symmetric 
and in each case the Hamiltonian was constructed so that the position operator x or 
the field operator (p is always multiplied by i. 

An experimental signal of a complex Hamiltonian might be found in the context 
of condensed matter physics. Consider the complex crystal lattice whose potential is 
V{x) = isin X. While the Hamiltonian H = p^ + isin x is not Hermitian, it is VT- 
symmetric and all of its energy bands are real. However, at the edge of the bands the 
wave function of a particle in such a lattice is always bosonic (27r-periodic) and, unlike 
the case of ordinary crystal lattices, the wave function is never fermionic (47r-periodic) 
[in]. Direct observation of such a band structure would give unambiguous evidence of 
a PT-symmetric Hamiltonian. 

There are many opportunities for the use of non-Hermitian Hamiltonians in the 
study of quantum field theory. Many field theory models whose Hamiltonians are 
non-Hermitian and PT-symmetric have been studied: PT-symmetric electrodynamics 
is a particularly interesting theory because it is asymptotically free (unlike ordinary 
electrodynamics) and because the direction of the Casimir force is the negative of that 
in ordinary electrodynamics [211] . This theory is remarkable because it can determine its 
own coupling constant. Supersymmetric PT-symmetric quantum field theories have also 
been studied [21]. A scalar quantum field theory with a cubic self-interaction described 
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by the Lagrangian C = ivr^ + ^(Vip)'^ + |/i^V^^ + gy^^ is physically unacceptable because 
the energy spectrum is not bounded below. However, the cubic scalar quantum field 
theory that corresponds to if in with e = 1 is given by the Lagrangian density 

This is a new, physically acceptable quantum field theory. 

We have found that PT-symmetric quantum field theories exhibit surprising and 
new phenomena. For example, consider the theory that corresponds to ii in (Q) with 
e = 2, which is described by the Lagrangian density 

c = y + uvvy + '^f^v-\9v'. (90) 

For example, for g sufficiently small this —g^'^ theory possesses bound states (the 
conventional gip'^ theory does not because the potential is repulsive). The bound states 
occur for all dimensions < D < 3 [5?, but for purposes of illustration we describe 
the bound states in the context of one-dimensional quantum field theory (quantum 
mechanics). For the conventional anharmonic oscillator, which is described by the 
Hamiltonian 

H = ]^p' + lm'x^ + -^gx^ > 0), (91) 

the small-^f Rayleigh-Schrodinger perturbation series for the /cth energy level Ek is 



k + - + -{2^ + 2k + l)v + 0(z/2) 
2 4 / \ / 



0+), (92) 



where v = g/{4m^). The renormalized mass M is defined as the first excitation above 
the ground state: M = - ~ m[l + 3u + ©(z/^)] as z/ 0+. 

To determine if the two-particle state is bound, we examine the second excitation 
above the ground state using We define 

B2 = E2- Eor^m[2 + 9u + 0{u^)] (z/ ^ 0+). (93) 

If B2 < 2M, then a two-particle bound state exists and the (negative) binding energy is 
B2 — 2M. If B2 > 2M, then the second excitation above the vacuum is interpreted as an 
unbound two-particle state. We see from (|^ that in the small-coupling region, where 
perturbation theory is valid, the conventional anharmonic oscillator does not possess a 
bound state. Indeed, using WKB, variational methods, or numerical calculations, one 
can show that there is no two-particle bound state for any value of g > 0. Because 
there is no bound state the gx'^ interaction may be considered to represent a repulsive 
force. Note that in general, a repulsive force in a quantum field theory is represented 
by an energy dependence in which the energy of a two-particle state decreases with 
separation. The conventional anharmonic oscillator Hamiltonian corresponds to a field 
theory in one space-time dimension, where there cannot be any spatial dependence. In 
this case the repulsive nature of the force is understood to mean that the energy B2 
needed to create two particles at a given time is more than twice the energy M needed 
to create one particle. 
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The perturbation series for the non-Hermitian, T'T-symmetric Hamiltonian 

H = ^p^ + ^m''x''-^gx^ (g > 0), (94) 

is obtained from the perturbation scries for the conventional anharmonic osciUator by 
replacing p — > —i>. Thus, while the conventional anharmonic oscillator does not possess 
a two-particle bound state, the T'T-symmetric oscillator does indeed possess such a 
state. We measure the binding energy of this state in units of the renormalized mass 
M and we define the dimensionless binding energy A2 by 

A2 = {B2 - 2M)/M ~ -3u + 0{u^) {v 0+). (95) 

This bound state disappears when v increases beyond v = 0.0465. As u continues to 
increase, A2 reaches a maximum value of 0.427 ai f — 0.13 and then approaches the 
limiting value 0.28 as — > 00. 

In the T'T-symmetric anharmonic oscillator, there are not only two-particle bound 
states for small coupling constant but also /c-particle bound states for ell k > 2. The 
dimensionless binding energies are 

Ak= (Bk-kM)/M -3k(k-l)u/2 + 0(u'^) (u^0+). (96) 

The coefficient of u is negative. Since the dimensionless binding energy becomes negative 
as 1/ increases from 0, there is a /c-particle bound state. The higher multiparticle bound 
states cease to be bound for smaller values of u; starting with the three-particle bound 
state, the binding energy of these states becomes positive as 1/ increases past 0.039, 
0.034, 0.030, and 0.027. 

Thus, for any value of u there are always a finite number of bound states and 
an infinite number of unbound states. The number of bound states decreases with 
increasing u until there are no bound states at all. There is a range of i> for which there 
are only two- and three-particle bound states, just like the physical world in which one 
observes only states of two and three bound quarks. In this range of u if one has an 
initial state containing a number of particles (renormalized masses) , these particles will 
clump together into bound states, releasing energy in the process. Depending on the 
value of v, the final state will consist either of two- or of three-particle bound states, 
whichever is energetically favored. There is a special value of u for which two- and 
three-particle bound states can exist in thermodynamic equilibrium. 

How does a g(p^ theory compare with a g(p'^ theory? A g(p^ theory has an attractive 
force. Bound states arising as a consequence of this force can be found by using the 
Bethe-Salpeter equation. However, the gcp^ field theory is unacceptable because the 
spectrum is not bounded below. If we replace g by ig, the spectrum becomes real and 
positive, but now the force becomes repulsive and there are no bound states. The same 
is true for a two-scalar theory with interaction of the form ig'-fP'Xi which is an acceptable 
model of scalar electrodynamics that has no analog of positronium. 

Another feature of PT-symmetric quantum field theory that distinguishes it from 
conventional quantum field theory is the commutation relation between the V and C 
operators. If we write C = + iCi, where Cr and Ci are real, then CrP = PCr and 
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CiV = —VCi. These commutation and anticommutation relations suggest the possibihty 
of interpreting PT-symmetric quantum field theory as describing both bosonic and 
fermionic degrees of freedom, an idea analogous to the supersymmetric quantum 
theories. The distinction here, however, is that the supersymmetry can be broken; 
that is, bosonic and fermionic counterparts can have different masses without breaking 
the VT symmetry. Therefore, another possible observable experimental consequence 
might be the breaking of the supersymmetry. 

8. Concluding Remarks 

We have described an alternative to the axiom of standard quantum mechanics that 
the Hamiltonian must be Hermitian. We have shown that Hermiticity may be replaced 
by the more physical condition of VT (space-time reflection) symmetry. Space-time 
reflection symmetry is distinct from the conventional Dirac condition of Hermiticity, so 
it is possible to consider new quantum theories, such as quantum fleld theories whose 
self- interaction potentials are ig(p^ or —g^p^. Such theories have previously been thought 
to be mathematically and physically unacceptable because the spectrum might not be 
real and because the time evolution might not be unitary. 

These new kinds of theories are extensions of ordinary quantum mechanics into 
the complex plane; that is, continuations of real symmetric Hamiltonians to complex 
Hamiltonians. The idea of analytically continuing a Hamiltonian was flrst discussed by 
Dyson, who argued heuristically that perturbation theory for quantum electrodynamics 
diverges [2,3J. Dyson's argument involves rotating the electric charge e into the complex 
plane e — > ie. Applied to the anharmonic oscillator (jHU, Dyson's argument goes as 
follows: If the coupling constant g is continued in the complex-g plane to —g, then 
the potential is no longer bounded below, so the resulting theory has no ground state. 
Thus, the ground-state energy Eo{g) has an abrupt transition at g = 0. As a series 
in powers of g, EQ{g) must have a zero radius of convergence because Ef){g) is singular 
ai g = 0. Hence, the perturbation series must diverge for all g 0. The perturbation 
series does indeed diverge, but this heuristic argument is flawed because the spectrum of 
the Hamiltonian (jMj) that is obtained remains ambiguous until the boundary conditions 
that the wave functions must satisfy are specifled. The spectrum depends crucially on 
how this Hamiltonian with a negative coupling constant is obtained. 

There are two ways to obtain H in ()94|) . First, one can substitute g = \g\e^^ into 
()91|1 and rotate from = to = tt. Under this rotation, the ground-state energy EQ{g) 
becomes complex. Evidently, E^i^g) is real and positive when (7 > and complex when 
g < Note that rotating from ^ = to 6' = — vr, we obtain the same Hamiltonian as 
in but the spectrum is the complex conjugate of the spectrum obtained when we 
rotate from ^ = to ^ = vr. Second, one can obtain ()94|) as a limit of the Hamiltonian 

H = \v' + \m^x^ + \gx\ixy {g > 0) (97) 
as e : 2. The spectrum of this Hamiltonian is real, positive, and discrete. The 
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spectrum of the limiting Hamiltonian (j94p obtained in this manner is similar in structure 
to that of the Hamiltonian in (j^?T|) . 

How can the Hamiltonian ()94j) possess two such astonishingly different spectra? 
The answer lies in the boundary conditions satisfied by the wave functions 0„(x). In 
the first case, in which 6 = sugg is rotated in the complex-t; plane from to vr, tpn{x) 
vanishes in the complex-x plane as |x| —>■ oo inside the wedges — vr/S < argx < and 
—Att/3 < aigx < —71. In the second case, in which the exponent e ranges from to 2, 
4>n{x) vanishes in the complex-x plane as |x| ^ oo inside the wedges — vr/S < argx < 
and — TT < argx < — 27r/3. In this second case the boundary conditions hold in wedges 
that are symmetric with respect to the imaginary axis; these boundary conditions enforce 
the VT symmetry of H and are responsible for the reality of the energy spectrum. 

Apart from the spectra, there is yet another striking difference between the two 
theories corresponding to H in ()94|1 . The one-point Green's function Gi{g) is defined 
as the expectation value of the operator x in the ground-state wave function (t)o{x), 



where C is a contour that lies in the asymptotic wedges described above. The value 
of Gi{g) for H in (jMj) depends on the limiting process by which we obtain H. If we 
substitute g = g^e^^ into the Hamiltonian (j^Tj) and rotate from ^ = to = tt, we find 
by an elementary symmetry argument that Gi{g) = for all g on the semicircle in the 
complex-(yf plane. Thus, this rotation in the complex-g plane preserves parity symmetry 
{x — > —x). However, if we define H in ()94|) by using the Hamiltonian in ()97p and by 
allowing e to range from to 2, we find that Gi{g) ^ 0. Indeed, Gi{g) ^ for all 
values of e > 0. Thus, in this theory VT symmetry (reflection about the imaginary 
axis, X — i> —X*) is preserved, but parity symmetry is permanently broken. 

Finally, we point out that the "wrong-sign" field theory described by the Lagrangian 
density (j90|l is remarkable because, in addition to the energy spectrum being real and 
positive, the one-point Green's function (the vacuum expectation value of the field is 
nonzero 24] . Furthermore, the field theory is renormalizable, and in four dimensions is 
asymptotically free (and thus nontrivial) [22] • Based on these features, we believe that 
a —gf^ quantum field theory the theory may provide a useful setting to describe the 
dynamics of the Higgs sector in the standard model. 
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